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Abstract 

Recently, diagrammatic extensions of dynamical mean field theory (DMFT) 
have been proposed for including short- and long-range correlations beyond 
DMFT on an equal footing. We employ one of these, the dynamical vertex ap¬ 
proximation (DFA), and study the two-dimensional Hubbard model on a square 
lattice. We define two transition lines in the phase diagram which correspond, 
respectively, to the opening of the gap in the nodal direction and throughout the 
Fermi surface. Our self-energy data show that the evolution between the two 
regimes occurs in a gradual way (crossover) and also that at low enough tem¬ 
peratures the whole Fermi surface is always gapped. Furthermore, we present 
a comparison of our DFA calculations at a parameter set where data obtained 
by other techniques are available. 
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1. Introduction 

Dynamical mean field theory (DMFT) has been a big step forward 

for the calculation and understanding of strongly correlated electron systems. 
It includes a major part of the electronic correlations: the local ones. In the 
vicinity of phase transitions or in one- or two dimensions non-local correlations 
are however essential. To include non-local correlations but to keep the local 
DMFT correlations at the same time, cluster extensions of DMFT have been 
developed such as the dynamical cluster approximation (DCA) and cluster 
DMFT [3, 0 . However numerical restrictions regarding the size of the clusters 
only allow to treat short range correlations this way. 

Hence, as an alternative and to deal with short- and long-range correlations 
on the same footing, diagrammatic vertex extensions of DMFT have been pro¬ 
posed more recently. The first such approach has been the dynamical vertex 
approximation (DFA) [^, followed by the dual fermion (DF) approach Q, the 
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one-particle irreducible approach (IPI) [1^, and the merger of DMFT with the 
functional renormalization group (DMF^RG) [HI- All of these approaches start 
with a local two-particle vertex [ij, [3 calculate from this non-local corre¬ 
lations beyond DMFT. This way, non-local correlations and associated physics 
are obtained, both, on the two-particle level as e.g. for the susceptibilities and 
on the one-particle-level as e.g. for the momentum-dependence of the self en¬ 
ergy. The difference between the various approaches lies in which two-particle 
vertex is taken: the fully irreducible one (full DTA), the irreducible one in a cer¬ 
tain channel (ladder DTA), the one-particle irreducible (IPI and DMF^RG) or 
the reducible one (DF). Then, depending on the approach, Feynman diagrams 
are constructed from the full Green function G or from the difference between 
G and the local Green function Gioc; and different kind of diagrams are taken: 
parquet, ladder, 2nd order or the ones generated by an RG flow. 

Hitherto, the diagrammatic extensions of DMFT have been applied to simple 
models, in particular the one-band Hubbard model, though the concept of ab 
initio calculations with DFA has also been proposed [1^. Physical highlights so- 
far have been the calculation of the critical exponents of the three-dimensional 

, quantum criticality 
. It was also possible 
to show that the paramagnetic phase of the s quar e lattice Hubbard model is 
always insulating at low enough temperatures T 23| , i.e. that the whole metallic 
side of the Mott-Hubbard metal-to-insulator transition as described by DMFT 
is completely washed out by extended long-range spin fluctuations in 2D. 

In this paper we recapitulate the DFA method in Section[2] for further details 
the reader is referred to [l^ and . Results for the two-dimensional Hubbard 
model are presented in Section [H Beyond [1^, we analyze further aspects 
of the transition from the high-T paramagnetic metallic phase to the low-T 
paramagnetic insulator. In particular, we start by identifying two crossover 
lines in the phase diagram: first the self-energy at the antinodal k = (tt, 0) 
turns insulating; at a lower T also nodal momentum k = (7r/2,7r/2) shows an 
insulating behavior. Specifically, our self-energy data demonstrates that the 
evolution between the two regimes takes place gradually by decreasing T for 
all the U values analyzed. In Section 01 we also present, as a benchmark, the 
comparison of DFA with DF and DGA. Section [5] summarizes our results. 


Hubbard model [Hill and the Falicov-Kimball model 
[l9l | , pseudogaps physics [1^ 211, and superconductivity 


2. Dynamical vertex approximation 

The idea of DFA is a resummation of Feynman diagrams in terms of the 
locality of the diagrams, not in orders of U. The first, one-particle level of 
this resummation is DMFT, which approximates the self energy, which is the 
one-particle fully irreducible vertex, to be local. Here, irreducible means that 
cutting one Green function line does not separate any self energy diagram into 
two pieces. Such reducible contributions are then generated by the Dyson Eq.: 

G(j.,k) = [l/Go(i^,k)-E(j.,k)]-i, (1) 
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Figure 1: Graphical representation of the Dyson equation. 


connecting the Green function G(j^, k) at Matsubara frequency v and momentum 
k (double blue line in Fig.[T|), its non-interacting counterpart Go(j^, k) = 

^ — Ck — k)] (single blue line) and the self energy I](i/, k) (red circles). 

The next, two-particle level then assumes the two-particle fully irreducible 
vertex to be local. This is DTA as employed nowadays. Before we discuss this 
approach in further detail, let us mention that, in principle, the fully irreducible 
n-particle vertex can be assumed to be local. This way more and more contri¬ 
butions are generated, and for n ^ oo the exact solution is recovered. However 
with increasing n the numerical effort also explodes, and going beyond n = 2, 
or possibly n = 3, does not seem to be feasible. 

Let us hence now focus on the two-particle level. The full two-particle vertex 
F consists of all connected diagrams. Some of these diagrams are two-particle 
fully irreducible, i.e., cutting two Green function lines G does not separate the 
diagram into two pieces. We denote these fully irreducible diagrams by A. Other 
(reducible) diagrams separate into two pieces. (Here and in the following we 
consider skeleton diagrams, i.e., Feynman diagrams in terms of the interacting 
Green function G instead of those in terms of Gq.) 

The reducible diagrams of F can be further classified. If we cut two G 
lines, there are three ways how F separates into two pieces: particle-hole ($p?i), 
particle-hole transversal ($^) and particle-particle (5*pp) reducible diagrams, 
see Fig. [2j One can prove, using particle-conservation, that each diagram is 
actually either fully irreducible or reducible in exactly one channel, so that 

T(1234) = A(1234) + %h{12M) + $^(1234) + 4>pp(1234). (2) 

Here, we have introduced a short-hand notation: 1 represents a momentum- 
frequency-spin coordinate 1 = (k, i/, cr), 2 = (k -|- q, z/ -|- w, cr), 3 = (k' -f q, -f 
w, cr'), 4 = (k'j z/', a'). Eq. (I2|) is known as the parquet equation. 

Eq. ([2|) has, for a given A, four unknowns: F and three <I>j.’s r G {ph,ph,pp}. 
Next, we hence first define the irreducible diagrams in a certain channel r: 

r^(1234) = F(1234)-$r(1234) . (3) 

From Fr, we can then calculate the reducible diagrams by repeatedly connecting 
them by two Green functions, i.e., formally <I>r = FrGGFr-l-rrGGrrGGrr-l-.... 
Inserting this as a recursion into Eq. ©, yields the Bethe-Salpeter Eqs. in the 
three channels (using Einstein’s summation convention), cf. Fig. |6l 

F(1234) = rp/,(1234)-FF(I22'T)G(3'2')G(l'4')rph(4'3'34) (4) 

= r^(1234)-f F(2'233')G(2'T)G(3'4')r^(lT4'4) (5) 

= rpp(1234)-FF(4'22'4)G(2'3')G(l'4')rpp(13'3l'). (6) 
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Figure 2: Graphical representation of the parquet equation: the full vertex F contains fully 
two-particle irreducible diagrams A and two-particle reducible diagrams in the particle-hole 
^phi transversal particle-hole and particle-particle reducible channel ^pp. 



Figure 3: Bethe-Salpeter Eq. 10 in the particle-hole (ph) channel. 


If we now substitute rr(1234) = f (1234) — (1234) in the three Bethe-Salpeter 
equations, we have four equations [Eqs. (O, (|1]), (IS|), dH])], which can be resolved 
for the four unknowns [F, ^ph, 4*^, and d^pp]. 

There is however one further complication, namely the internal Green func¬ 
tion lines G in the Bethe-Salpeter Eqs. These G’s can be calculated from E via 
the Dyson Eq. O; and E in turn can be calculated exactly from F through the 
Heisenberg Eq. of motion (also called Schwinger-Dyson Eq. in this context): 

E(14) = -t/(12'3'l')G(l'4')G(23')G(2'3)F(4'234) 

-f[/(1234)G(23) - C/(1432)G(23) (7) 


That is the solution of the exact parquet Eqs. involves six equations [(E®, 
®, ®, 0, ®] and six unknowns [F, 4)p/i, 4)^, 4)pp, E, G]. If the exact A 
were known, this would allow us to calculate the exact vertex F, E, and G. If 
not, we still can do approximations for A. The so-called parquet approximation 
assumes A = U 2^; in DTA we consider all Feynman diagrams for A in all 


orders of C/, but only take their local contribution. Let us note that the locality 
assumption for A is much better fulfilled than that for E. Even for the two 
dimensional model considered in this paper. DCA calculations have shown A 
to be essentially k-independent while E and F are strongly k-dependent. 

So-far we have discussed the full (or parquet) DTA. Often (also for the half- 
hlled Hubbard model discussed in this paper) spin-fluctuations dominate over 
charge fluctuations at low energies. In this case, the particle-hole and particle- 
hole transversal channels dominate over the particle-particle channel. Hence, 
we can neglect 4)pp. Then, th e 4>n fr and 4)^ channels decouple, and are even 
related by crossing-symmetry [l^, so that the calculation of the two Bethe- 
Salpeter Eqs. dH), ([6|) (or even one because of the crossing symmetry with a 
local Tp/j as a starting point) is sufficient. This ladder-DTA is numerically much 
easier and the two channels modulo the double counting can be combined to 
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yield the full vertex in Eq. ([2|) 13. 2^. 


3. Metal-insulator crossover in the 2D Hubbard model 


Let us now turn to the ladder-DFA results for the half-filled Hubbard model 
on a square lattice with nearest neighbor hopping t and local Coulomb repulsion 
U, given by the Hamiltonian 


H = -t 


(d'>o- 


^ia^ja 






( 8 ) 


Here, (Co-) creates (annihilates) an electron on lattice site i with spin a, (ij) 
sums each nearest neighbor pair once, and energies are from now on measured 
in units of the half-bandwidth D = 4t. 


The main physical result of 23| is that the low-temperature paramagnetic 
phase (above the T = 0 antiferromagnet) is always insulating for all U > 0. In 
other words, the metal-insulator transition is at £4 = 0 instead of a finite value 
Uc > 0 which was previously concluded from cluster extensions of DMFT, see, 
C-g-j 2M- As a function of temperature T however, there is a crossover from 
the low-T paramagnetic insulator to a high-T paramagnetic metal at small U’s. 
Here, we analyze in greater detail how this crossover occurs. 

Fig. |T]shows the imaginary part of the self-energy for the nodal and antinodal 
points of the Fermi surface. For U = 0.5 and high temperatures, T > 0.017, the 
self-energy shows clearly metallic behavior at every point on the Fermi surface 
from the nodal to the antinodal one. Upon cooling along the gray arrow in Fig. 
01 first the antinodal point of the Fermi surface [k = (tt, 0)] shows a downturn 
for 0 in Fig. 01 i.e., an insulating behavior below T = 0.017. 

For 0.017 > T > 0.0125 however, the nodal point of the Fermi surface 
[k = {-^12, 7r/2)] does not show this insulating behavior yet. Only for T < 0.0125 
this and all other points of the Fermi surface show an insulating self energy. In 
other words, for T < 0.0125 we have an insulator with the whole Fermi surface 
gapped, whereas for 0.017 > T > 0.0125 we have a pseudogap with only parts 
of the Fermi surface gapped. 

Hence, upon cooling along the grey arrow in the phase diagram Fig. 01 
we hence first cross the dashed red line which marks the temperature where 
the first point of the Fermi surface [k = (tt, 0)] shows insulating behavior. At 
lower temperatures, we cross a second (solid-red) line at which the whole Fermi 
surface gets gapped. The red-shaded region inbetween hence has a pseudogap. 

The analogous analysis for U = 1.0 can be found in the lower panels of Fig. 
01 It corresponds to the second gray arrow in Fig. 01 In this case, we identify a 
larger pseudogap region at higher temperatures: for 0.05 > T > 0.025 the nodal 
self energy is still metallic whereas the antinodal one is already insulating. 


4. Comparison of numerical techniques for the 2D Hubbard model 

As mentioned in the Introduction, for the treatment of the two-dimensional 
Hubbard model there exists a variety of numerical techniques, that aim at 


5 






0 

- 0.01 
- 0.02 

d -0.03 

w 

E -0.04 
-0.05 
-0.06 

0 0.1 0.2 0.3 0.4 0.5 0 0.1 0.2 0.3 0.4 0.5 




0 

- 1 - 1 - 1 -1 - 1 -1-1- 

0 

- 1 - 1 - 1-1 - 1 - 1 -1— 

-0.1 


-0.1 


-0.2 

U=1.0 

7 k=(7i,0) 

-0.2 

7 W=(kI2,tiI2) ■ 

-0.3 

7 T=0.05 — T=0.017 T=0.01 —' 

/ T=0.029 — T=0.0125 — DMFT 

-0.3 

■J 

-0.4 

1 I I 1 I 1 1 

-0.4 

■ 1 1 1 I I 1 


0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 


Figure 4: Imaginary parts of the DFA self-energy S vs Matsubara frequency Un for the half- 
filled Hubbard model at the antinodal (left) and nodal (right) point of the Fermi surface at 
U = 0.5 (upper panels) and U = 1.0 (lower panels) and different temperatures. The DMFT 
results at T = 0.01 is provided for comparison. 


including correlations beyond the purely temporal but local ones covered by 
DMFT. Since there is no exact solution of this model in 2D, benchmarks be¬ 
tween different methods are of fundamental importance. Fig. |6] shows a self¬ 
energy comparison of our DMFT and ladder-DFA data with the DCA and DF 
results extracted from the extensive numerical review [ 2 ^ by J. LeBlanc et al. 
for U = 2.0 and T = 0.125. These parameters correspond to the black triangle 
in the upper right corner of Fig. [S] 

In the main panel of Fig. |6] one can see the data for the antinodal point 
k = (tTjO). Because of the high temperature and large interaction all data sets 
display an insulating behavior, also DMFT. However, the inclusion of non-local 
correlations (DCA, DF, DFA) further enhances the insulating tendencies; the 
self-energy downturn is more pronounced. The magnitude of this effect is similar 
in DF and DFA and somewhat larger in DCA. The overall trend is similar for 
the nodal point k = (7r/2,7r/2) shown as the inset of Fig. |6l In contrast to the 
weak coupling region which we focused on in Section[3l for this larger [/-value we 
have a Mott-Hubbard insulating behavior for all k points of the Fermi surface. 
This reflects the increased importance of local correlations in this regime; even 
DMFT yields resonable results. 
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Figure 5: DFA phase diagram for the half-filled Hubbard model on a square lattice. Due to 
strong antiferromagnetic fluctuations, we have a paramagnetic insulator for all C/ > 0 below 
the solid-red line 1^ : upon cooling (e.g. along the vertical arrows that mark the temperatures 
presented in Fig. 0 first the antinodal point [fc = (vr, 0)] turns insulating at the dashed-red 
line where the other k points still show a metallic self energy. That is, in the red-shaded 
region, we have a crossover with a pseudogap. The dashed-gray and solid-gray line indicates 
the Neel temperature Tjy in DMFT and DFA for 3D, respectively (for 2D Tjy = 0 in DFA). 
The triangle marks the parameters for the comparison in Fig. |6] 


5. Conclusion 


Due to long-ranged antiferromagnetic correlations, the paramagnetic phase 
of the two-dimensional Hubbard model is insulating at low enough temperatures, 
half-filling, and a lattice with perfect nesting such as the square lattice with 


nearest-neighbor hopping — even at an arbitrary small coupling U 23|. For 


small U and higher temperatures, on the other hand, we have a paramagnetic 
metal. In this paper, we analyzed the metal-to-insulator crossover in detail. We 
find that upon cooling, first a gap opens at the antinodal point before, at a lower 
temperature, the full Fermi surface is gapped. That is, there is a pseudogap. 
We have also compared our DFA results with DF and QMC where these are 
available, i.e., at higher temperatures and intermediate coupling strength. In 
this parameter regime, DMFT already gives a good description, and all of the 
three beyond-DMFT approaches show similar corrections to DMFT: non-local 
correlations systematically increase the insulating tendencies. 
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Figure 6: Comparison of the imaginary parts of DMFT iblue squares), DFA (red dots), DCA 
(grey triangles from 1281 ) and DF (grey stars from l28l ) self-energies for the half-filled 2D 
Hubbard model with U = 2.0 and T = 0.125. The main plot shows the antinodal point 
(fc = ('7r,0)) and the inset the nodal point (fc = (^72, 7 r/ 2 )) of the Fermi surface. For the 
latter, no DF/DCA data are available for comparison. 
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